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INVESTIGATION OF WO BLUFF SHAPES IN AXIAL FREE FLIGHT 
OVER A MACH NUMBER RANGE FROM 0.35 TO 2.15 
By Lucille C. Coltrane 


SUMMARY 


A fineness -ratio-2 . 71 right circular cylinder and a fineness -ratio- 
4.0 circular cylinder with a rounded nose and a flared afterbody have 
been tested in free flight over a Mach number range of 0.35 to 2.15 and 

a Reynolds number range of 1 X 10^ to 12 x 10^. Time histories, cross 
plots of force coefficients, rolling velocity, and longitudinal -force 
coefficient are presented for both cylinders. In addition, cross plots 
of moment coefficients and plots of the normal -force curve slope and 
the aerodynamic center are presented for the fineness -ratio- 2.71 cylinder. 
The average aerodynamic center of the right circular cylinder moved rear- 
ward with decreasing speeds until at the subsonic Mach numbers it remained 
approximately constant and comparisons of the drag data of this test with 
wind tunnel and other free -flight data show good agreement. An appreciable 
decrease in drag was observed when the data of the present test of the 
rounded nose cylinder were compared with data of a right circular cylinder 
of a similar configuration. 


INTRODUCTION 


The dynamic stability characteristics of shapes having very high 
drag and very low lift are of interest to designers of various aerodyna- 
mic stores (see refs. 1 and 2) and shapes such as a reentry missile. One 
configuration which meets these requirements is the circular cylinder in 
axial flow. An investigation is being conducted by the Langley Pilotless 
Aircraft Research Division with the rocket-boosted free -flight -model 
technique to obtain stability data on bluff shapes in axial free flight. 

The present investigation was conducted to extend the data in ref- 
erence 3 to higher Mach numbers and to determine the effect of a change 
in nose shape on the aerodynamic characteristics of a bluff shape . For 
this purpose flight tests were made of two additional cylinders, a right 
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circular cylinder with a fineness ratio of 2.71 of the same design as 
the fineness -ratio - 2. 56 cylinder of reference 3 ant ^ a circular cylinder 
with a rounded nose and a flared afterbody of the same design except for 
nose shape as the fineness -ratio-4 .0 cylinder of reference 3 . The 
results of these flight tests are presented in this paper. These tests 
covered a Mach number range from 0.35 "to 2.15 and a Reynolds number 

range from 1 X 10^ to 12 x 10^ based on cylinder diameter. The free- 
flight tests were conducted at the Langley Pilotless Aircraft Research 
Station at Wallops Island, Va. 


SYMBOLS 


The data are presented relative to the body axis system and the 
positive directions of the force coefficients, moment coefficients, and 
angular velocities are shown in figure 1. The various symbols used 
throughout the paper are defined as follows: 


a accelerometer reading, g units 


longitudinal -force coefficient, 


a l,cg 


w/s 

q. 


m 


ly /•• • • N 

pitching -moment coefficient, (0 - \|rqp) 

qSd 


dC. 


"m, 


m 

n • 

2V 


'N 


normal -force coefficient, a n,cg 


w/s 


dC 


"N, 


N 


a da 


'n 


yawing -moment coefficient, — — (^ + 9<p) 

qSd 


lateral -force coefficient, a^ 
diameter of cylinder, ft 


w/s 

q 
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g 

I 

T X’ X Y’ 
k 

l 

M 

m 

q. 

R 

r 

S 

t 

V 

w 

x 


0 

A 0 
A A 

U 

P 

T 

V 


acceleration due to gravity, ft/sec 2 

O 

moment of inertia, slug-ft 

moments of inertia, slug-ft 2 


radius of gyration, 



ft 


length of cylinder, ft 
Mach number 
mass, slug 

dynamic pressure, lb/sq ft 
Reynolds number based on cylinder diameter 
radius of cylinder nose, in. 
cross-sectional area of cylinder, sq ft 
time, sec 

free -stream velocity, ft/sec 
weight of cylinder, lb 

distance along cylinder from nose, ft, or, when used with a 
subscript, distance from center of gravity, positive 
forward, ft 

angle of pitch 

nonrolling damping constant, l/sec 
damping constant due to roll, l/sec 

2i— 

relative -density factor, 

pSd 

air density, slugs/cu ft 
angle of roll, radians 
angle of yaw, radians 
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“0 


0 - \jfcp 


basic oscillation frequency, radians/sec 

component of total pitch frequency resulting directly from 
roll, radians/sec 

a 2 ” a n ]_ 

pitching acceleration, g — - — , radians/sec^ 

x n,2 " ^,1 


n,cg 


normal acceleration, 


^jl-^n^ ~ a n, 2 x n,l 
x n,2 - %,1 


, g units 


• • • • cL-t 2 _ 1 / p 

\|r + 9® yawing acceleration, g — - — , radians/sec 

x t,2 " x t,l 


t,cg 


transverse acceleration. 


^•t , l x t , 2 ~ a t , 2 x t , 1 
x t,2 " x t,l 


, g units 


Subscripts : 

ac aerodynamic center 

eg center of gravity 

l longitudinal 

n normal 

t transverse 

tot total 

1 forward end of cylinder 

2 rear end of cylinder 

A dot above a symbol indicates time rate of change of symbol, for 
example, 0 = 80/8 t. 


MODEIS 


The physical characteristics of the models are presented in table I, 
and drawings of the models are shown in figure 2. Photographs of the 
models and model booster arrangements are presented in figure 3. 
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The fineness-ratio-2.71 model (designated model I) was a right cir- 
cular cylinder with a diameter of 8 inches and contained two small rockets 
which were mounted normal to the longitudinal axis and ahead of the center 
of gravity to give a yaw disturbance. The center of gravity of this 
model was located at 34.3 percent of the cylinder length behind the nose. 
The fineness-ratio-4.0 model (designated model II) was a circular cylinder 
with a diameter of 8.0 inches, a 13 ° flared afterbody, and a rounded nose 
and contained no small rockets . The afterbody and nose lengths are shown 
in figure 2 and the contour ordinates for the nose are shown in table II. 
The center of gravity of this model was located at 24.5 percent of the 
model length behind the nose. 

The cylindrical portions of the models were constructed of steel and 
were covered with a fiber-glass plastic shell. The forward portion of 
model I and the nose of model II were machined from solid steel . The 
afterbody flare of model II was constructed of laminated wood. 


INSTRUMENTATION 


Model instrumentation consisted of an NACA six-channel telemeter 
which transmitted data from six accelerometers located as follows : one 

normal and one transverse accelerometer in the forward end of the 
cylindrical portion of the model, one normal and one transverse accel- 
erometer in the rear of the cylindrical portion of the model, and two 
longitudinal accelerometers, one for high range and one for low range, 
behind the center of gravity. A measure of the signal strength trans- 
mitted from the loop antenna provided an indication of the roll rate of 
the cylinders since the strength of the signal varied with the cylinder 
roll position. 

Ground instrumentation included a CW Doppler radar unit to measure 
the velocity of the model, a modified SCR 5^4 tracking radar set to 
determine the flight path and a rollsonde receiver used as an additional 
measure of the rolling velocity. Fixed and tracking motion -picture 
cameras were used to observe the model during the first portion of the 
flight. Atmospheric data were obtained from a rawinsonde released immedi- 
ately before model flight. 


TESTS AND ANALYSIS 


Tests 


The models were ground launched at an angle of 70° from the 
zontal by means of a mobile launcher as shown in figure 3(c). A 
fuel Cajun rocket motor boo s^gd the models^Jb o maximum velocity. 
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flaps were incorporated into the booster to increase the separating force 
between the model and booster at booster burnout. For model I, the drag 
flaps operated as designed, but, for model II, motion pictures showed the 
drag flaps opened prematurely, came off, and damaged the booster fins. 

The damaged fins probably increased the drag of the booster and, although 
separation was delayed about 0.5 second, a clean separation was experi- 
enced. In addition, the damaged booster fins caused the model-booster 
combination to experience a high rolling velocity immediately before sepa- 
ration. This high rolling velocity persisted throughout the flight of the 
model. Tracking radar showed that the models followed an approximately 
parabolic flight path. 


Analysis 

As was done in reference A, the method of analysis used herein to 
determine the pitching -moment coefficients consists of using the data 
obtained from normal accelerometers located at two positions in the model. 
When the normal-force equation is combined with the equation describing 
the normal -accelerometer reading at any point along the x-axis of the 
model, it can be shown that 


(=n,2 - =n,l) = e - W (1) 

*n,2 ' *n,l v ' 

The pitching -moment equation may he written as 

qSdC mtQt = lye - (I Z - I X )<tf (2) 

if ~ 0. Now, if Iy = l£ as is the case for the models reported 

herein (see table I), 


iSdCm^^ = Iy(9 - qxjf ) + I x q^ 

and then, if Iycp^ is small with respect to Iy(0 - <pjf). 


(3) 


C ">tot - - «> 


(4) 


For model I, equation (4) may be used directly since the roll rate 
was found to be small and the term IyCpj/' cgn be neglected. 


° 
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or 


For the aerodynamic -center location, where 


_ dC m d x cg 
ac dC N l l 


x = ££n £l + 
ac d.Cy l + l 


the values for 


E 

dC- 


dC r 


and were taken directly from the polars of 


■m 


against Cjvj and C n against Cy at or Cy near zero. 


In order to determine the normal -force coefficient curve slope as 
presented herein for model I and for the fineness -ratio-2. 56 cylinder 
of reference 3 , equation ( 42 ) in reference 5> which states that 


'N. 


-2nhr - -r 


_d_ 

2 V 


a 


if! 

«, 2 + V 3 - M 2 + *» 2 1 A — ill 

- 1 2 


( 5 ) 


was simplified by omitting the damping terms. Analysis of the data that 
were obtained after the firing of the two pulse rockets showed that the 
last three terms within the brackets were small when compared with the 
term of the basic oscillation frequency co(^. (See table III.) There- 
fore, they were neglected throughout the calculation. Then by the 
relation 


s = 


xg 

d 


(k/d) 


I 1 

d ) 

dC N 


dC„ 


( 6 ) 


dC„ 


it is possible to express Cjg- in terms of as follows 


a 


m 
dC N 


c, 
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(7) 


Tj-10 slopes — — and — — were equal for "the polars of model X . Since 
dCjj dCy 

variations of C N and a were probably nonlinear, this average C N ^ 

may be greater or less than values near Cjj = 0, depending upon the 
range of C^j encountered in flight. 


An indication of the dynamic stability following the disturbance 
given the model by the pulse rockets was obtained through the use of 
equation (43) in reference 5, where 



= 2 


!) 2 (¥ 


A 0 + 



(8) 


ACCURACY 


The instrument inaccuracies , estimated to be +2 percent of their 
calibrated range, are stated in coefficient form for representative Mach 
numbers as follows : 


Coefficient 

Fineness -ratio-2. 71 
model (model I) at 
Mach numbers of - 

Fineness -ratio -b . 00 
model (model II) at 
Mach numbers of - 

1.88 

1.26 

0.80 

2.00 

1.38 

0.80 

C N 

±0.039 

±0.094 

±0.250 

±0.038 

±0.085 

±0.309 

Cy 

±0.040 

±0.095 

±0.254 

±0.039 

±0.087 

±0.316 

c c 

±0.051 

±0.030 

±0.079 

±0.052 

±0.030 

±0.104 
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PRESENTATION OF RESULTS 


The variation of test Reynolds number, based on cylinder diameter, 
with test Mach number is presented in figure 4. The model flight paths 
are presented as plots of altitude against horizontal distance in fig- 
ure 5 and the variation of velocity and dynamic pressure with time is 
shown in figure 6. 

Time histories of the normal -force coefficient, lateral -force coef- 
ficient, and Mach number are presented in figures 7 and 8. Basic data 
cross plots of force and moment coefficients for model I are shown in 
figure 9* Examples of the variation of the normal -force coefficient with 
lateral-force coefficient for model II are shown in figure 10. An indi- 
cation of the rolling velocity against Mach number is shown in figure 11. 
For model I, as a function of Mach number, the variation of the average 
aerodynamic center in cylinder lengths behind the nose is presented in 
figure 12, and the normal -force curve slope is presented in figure 13 . 

The measured longitudinal -force coefficient for both models is plotted 
against Mach number and is presented in figure 14. 


A time history of M, Cjj, and Cy for model I is shown in fig- 
ure 7* A disturbance at separation (t = 2.85) caused the instruments to 
exceed their calibrated ranges at each oscillation peak for about 1 sec- 
ond. Throughout the flight a coupled motion with respect to the body axis 
system was ejqperienced by the model as the oscillations damped from large 
amplitudes at supersonic speeds to low amplitudes at subsonic speeds . 
Figure 8 presents a time history of M, Cjj, and Cy for model II. 


Plots of Cjj against Cy presented in figures 9 and 10 for various 

Mach numbers indicate the model motion. When the trim center is estimated, 
the time history of the resultant force vector can be obtained and the 
angular displacement of adjacent peaks gives an indication of the rolling 
rate. This use of the Cjj and Cy cross plots is discussed in detail in 

reference 5- Also, as stated in reference 5 > whether the missile is above 
or below roll resonance is immediately determined, since the motion above 
resonance is characterized by inside loops, and the motion below roll 


DISCUSSION OF RESULTS 


Time History 


Basic Data Plots 
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resonance by outside loops. Therefore, it can be seen that model I was 
below roll resonance and model II was above roll resonance; in fact, by 
comparing the plot of Cjj against Cy of figure 10 with a similar plot 

in figure 5 of reference 5, it is apparent that model II rolled at a rate 
of about three times the pitch frequency. This ratio is nearly the same 
as the ratio of the measured roll rate of model II and the nonrolling 
pitch frequency of a nearly similar fineness ratio 4.0 cylinder of ref- 
erence 3. Dashed lines in figure 9 indicate extrapolated values where 
the range of the instruments was exceeded but these data points were con- 
sidered accurate enough for obtaining an estimate of the roll rate. This 
estimate is shown with measured averages from the rollsonde in figure 11 
and good agreement is noted. 

Figure 9 shows that the variation of moment coefficient with force 
coefficient for model I was nonlinear in both the pitch and yaw planes, 
the greatest stable slope generally occurring at the zero force coeffi- 
cient. An attempt was made to determine the stability characteristics for 
model II from the disturbance received at separation; however, the scatter 
of data for the plots of C m against C N and C n against C Y pro- 
hibited the determination of a slope dC m /dC N to be used in equation (6). 
Thus neither Cm nor the aerodynamic -center location could be deter- 

mined. The relative importance of the various terms within the brackets 
of equation (5) for a model rolling about 100 radians/sec can be noted in 
table III for model II. 


Aerodynamic Center 

The variation of the aerodynamic center with Mach number for model I 
is shown in figure 12. The average aerodynamic center moved rearward as 
the speed decreased through the supersonic and transonic speeds but 
remained at a constant position of about 0.64 cylinder lengths behind 
the nose through the subsonic range. For comparison, a dashed-line curve 
faired from data of the fineness -ratio -2. 56 cylinder of reference 3 is 
also presented in figure 12. The comparison is considered good especially 
when the accuracies in the basic data are taken into account. 


Lift 

The variation of the normal -force curve slope with Mach number for 
model I is shown in figure 13 . It is seen that Cjj decreased with 

decreasing speeds until at subsonic Mach numbers it remained approximately 
constant. Also shown in figure 13 is a curve of Cm from faired data of 

the fineness -ratio-2. 56 model of reference 3 which shows good agreement. 
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With Cjj^ known, was calculated to be -2 6 per radian at 

M = 1.03 and 22 per radian at M = O.7O. Since the term Aq is small 

(see table III) and in these flight tests is difficult to obtain with 
good accuracy, these values should be considered as good estimates and 
not absolute values of the dynamic stability at these Mach numbers . 


Drag 

The variation of longitudinal -force coefficient with Mach number 
for the models of the present tests is presented in figure lk. The 
shaded areas on the figures indicate values of Cq obtained for all Cjj 

and Cy values of the tests . Also presented are the drag data for both 
models when Cjj < 0.1 and Cy < 0.1. For model I, comparison of these 

data with the drag of the fineness-ratio-2.56 cylinder of reference 3 
when Cjj < 0.1 and Cy < 0.1 and the drag data of reference 6 for the 

body of revolution where r = 0.00 inch shows good agreement. Compari- 
son of the data of model II with the data of the fineness -ratio-k.O cylinder 
with flared afterbody of reference 3 when Cjy < 0.1 and Cy < 0.1 indi- 
cates that a considerable reduction in drag throughout the Mach number 
range was obtained by the nose modification of the present test. The drag 
data of model II and the drag data of reference 6 (where r = 0.50 inch) 
show reasonable agreement at subsonic speeds . 


CONCLUDING REMARKS 


A fineness -ratio-2. 71 right circular cylinder and a fineness -ratio- 
k.O cylinder with a rounded nose and a flared afterbody have been tested 
in rocket-boosted axial free flight over a Mach number range of 0.35 to 

2.15 and a Reynolds number range of 1 X 10^ to 12 x 10^. When given a 
disturbance at supersonic speeds, the response of the fineness-ratio-2.71 
cylinder was a large -amplitude oscillation which damped to a low -amplitude 
sustained oscillation at subsonic speeds. The model experienced a low roll 
rate throughout the test. This value was less than the roll rate required 
for resonance. The moment-force curves were nonlinear, with greatest 
stable slope generally occurring at the zero force coefficient. The aver- 
age aerodynamic center moved rearward and the normal -force curve slope 
decreased with decreasing speeds until at the subsonic Mach numbers they 
remained approximately constant. Comparisons of the drag data of this 
test with wind-tunnel and other free -flight data show good agreement. 

The f ine ness -rati o -k. 0 cylinder experienced a high rate of roll throughout 
the test. This value was considerably more than the roll rate required 


C 


CONFIDENTIAL 


* 


12 


OOSFIBENTXitt 


• •• 


MCA RM L58A16 


for resonance. An appreciable decrease in drag was observed when the 
data of the present test of a rounded nose cylinder were compared with 
data of a right circular cylinder of a similar configuration. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 30, 1957- 


REFERENCES 


1. Walchner, Otto, Weigand, Heinrich, and Kroeger, Hemann: Dynamic 

Stabilizer for the Finless Bluff Bomb. Tech. Note WCRR 54-53, 

Wright Air Dev. Center, U. S. Air Force, Apr. 19, 1954. 

2. Carter, Howard S., and Lee, John B.: Investigation of the Ejection 

Release of Several Dynamically Scaled Bluff Internal Stores at Mach 
Numbers of 0.8, 1.39) and 1.98. NACA RM L 56 H 28 , 1956. 

3. McFall, John C., Jr.: Dynamic Stability Investigation of Two Right 

Circular Cylinders in Axial Free Flight at Mach Numbers From 0.4 
to 1.7 - Fineness -Ratio-2 . 56 Cylinder and Fineness -Ratio -4.0 Cylinder 
With Flared Afterbody. NACA RM L 56 L 28 , 1957* 

4. Peck, Robert F., and Coltrane, Lucille C.: Longitudinal Characteris- 

tics at Transonic and Supersonic Speeds of a Rocket-Propelled Air- 
plane Model Having a 60° Delta Wing and a Low Swept Horizontal Tail. 
NACA RM L 55 F 27 , 1955 . 

5 . Nelson, Robert L. : The Motions of Rolling Symmetrical Missiles 

Referred to a Body -Axis System. NACA TN 3737) 1956. 

6 . Polhamus, Edward C.: Effect of Nose Shape on Subsonic Aerodynamic 

Characteristics of a Body of Revolution Having a Fineness -Ratio of 
10.94. NACA RM L57F25, 1957- 


>| CONFIDENTIAL 


MCA RM L58AI6 



TABLE I.- PHYSICAL CONSTANTS FOR CYLINDERS TESTED 
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Constant 

Model I 

Model II 

W, lb 

Iy, slug-ft 2 . . . 

l y, SlUg -ft 2 . . . 

l z , slug -ft 2 . . . 

X cg /l 

d, ft 

1, ft 

S, ft 2 

97-25 

0.175 

0.910 

0.910 

0.3^3 

O.667 

1.805 

0.348 

108.5 

0.238 

1.430 

1.430 

0.245 

0.667 

2.667 

0.348 


TABLE II . - NOSE ORDINATES FOR MODEL II 


x, in. 

y , in. 

0 

1.600 

.01 

2.198 

.10 

2.765 

.20 

2.964 

.30 

3.086 

.40 

3.177 

• 50 

3-249 

.60 

3-309 

.80 

3.406 

1.00 

3-483 

1.50 

3.627 

2.00 

3-733 

2.50 

3.817 

3.00 

3.887 

3.50 

3.947 

4.00 

4.000 
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TABLE III.- VALUES USED TO DETERMINE SLOPE OF 
NORMAL-FORCE COEFFICIENT 


Mach 

number 

u>o 

\> 

AX 

Au 

^ - $ 
(1 - sr 


Model I 

1.05 

22.5 

-1.33 

0 

2.75 

0.224 

0.70 

14.0 

+0.57 

0 

1-95 

0.224 


Model II 

2.07 

37.0 

-2.33 

1.13 

95*3 

0.190 
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Figure 1.- Axes system showing positive directions of force and moment 
coefficients and angular velocities. 
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(a) Model I: Fineness -ratio -2.71 cylinder. 


^ 7.8k w 


*- k.o+j 

1 




52.0 




(b) Model II: Fineness-ratio-4.00 cylinder with a rounded nose 

and flared afterbody. 

Figure 2.- Drawings of cylinders tested. All dimensions are in inches. 
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(a) Model I. L- 57-1733 



(b) Model II. L- 96363 .I 

Figure Photographs of cylinders tested. 
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(c) Model II on booster in launching position. 


L- 57-1356.1 


Figure 3*- Concluded. 
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0 .4 .8 1.2 1.6 2.0 2.4 


M 



0 4 8 12 16 20 24 x 10 5 

Horizontal distance, ft 


Figure 5.- Flight paths of cylinders tested. 
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Figure 6.- Velocity and dynamic pressure of cylinders tested. 
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Figure 7 .- Time history of fineness-ratio-2.71 cylinder. Model I. 
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Figure 7*- Concluded. 
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Figure 8.- Time history of fineness-ratio-4.00 cylinder with flared 
afterbody and rounded nose. Model II. 
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Time, sec 


Figure 8.- Concluded. 
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(a) M = 2.09 to M = 1.66. (b) M = 1.66 to M = I. 34 . 

Figure 9« _ Basic data cross plots of force and moment coefficients for 
the f ineness-ratio-2 . 71 cylinder. Model I. The time sequence is 
indicated by the symbols O, □ , and A . 
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-2.5 -2.0 - 1.5 _1 *° ~»5 0 .5 1.0 1.5 2.0 2.5 

c » 

(c) M = 1.35 to M = 1.18. 

Figure 9»- Continued. 
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(d) M = 1.09 to M 


0.84. 


Figure 9*- Continued. 
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Figure 9;- Continued. 
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(g) M = O.58 to M = 0.49. 


(h) M = 0.49 to 


M = 0.3c 


Figure 9.- Concluded. 


CONFIDENTIAL 


30 


NACA RM L58A16 


• • • * 
• • 

• • 


• • • 
• • 


CC 


CONFIDENTIAL 


i 






Cy C Y 


Figure 10.- Examples of force coefficients for the fineness -ratio-4. 00 
cylinder with flared afterbody and rounded nose. Model II. The 
time sequence is indicated by the symbols O , □ , and A . 
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Figure 11.- Variation of rolling velocity with Mach number. 
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Figure 12.- Variation of aerodynamic center with Mach number. 



Figure 13.- Variation of normal -force coefficient slope with Mach 

number . 
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O Fineness-ratio-2.71 cylinder of this test. 

(Cc where C N < 0.1 and Cy-CO.l) 

□ Fineness -rati o-2. 56 cylinder of reference 3 . 
(Cc where CflCO.l and CyCO.l) 

A C c of reference 6 (where r - 0 . 00 ) 



Figure lL.- Variation of longitudinal -force coefficient with Mach 
number. Shaded area indicates measured C c for all C N and 

C Y values of this test. 
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